Topics in Guided Wave Propagation 
Through Gyromagnetic Media 

Part III — Perturbation Theory and Miscellaneous Results 

By H. SUHL and L. R. WALKER 

Some problems, complete discussion of which would be extremely difficult, 
are treated approximately by means of perturbation theory. Among these 
are the partially filled cylindrical waveguide, and the problem of ?nultiple 
internal reflections in a sav^ple of finite length filling the cross section of a 
cylindrical guide. Propagation in a ferriie between parallel planes, mag- 
netized along the propagation direction is discussed by the methods described 
in Part I. The paper concludes with an addendum to Part I — a numer- 
ical study of field patterns of the TEn-limit and TM n-limit mode for 
various dc magnetic fields. 

INTRODUCTION 

Parts I and II of this paper were devoted to a number of specific propa- 
gation problems, whose solutions, though frequently quite complicated, 
could be discussed with a reasonably modest investment of effort. Un- 
fortunately, not all of these problems pertain to situations met with in 
actual gyromagnetic devices. Actual devices frequently employ struc- 
tures whose performance could be predicted only as the result of lengthy 
computing programs. For example, the microwave gyrator using Faraday 
rotation usually employs a ferrite sample whose cross-section only partly 
fills that of the cylindrical waveguide. Although it is easy to formulate 
the corresponding equation for the propagation constant, the classifica- 
tion and survey, let alone the computation of solutions, would be very 
difficult to carry out. 1 

Thus, one must often be content with approximate results, and the 
bulk of the present paper is devoted to perturbation methods. These 
take as starting point a situation whose propagation problem is essen- 
tially solved. The small change in propagation constant due to a slight 
change in the original state of the sj'^stem is then calculated. The small 
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change of state may be the weak magnetization of an originally un- 
magnetized specimen occupying a substantial part of the structure, or 
the introduction of a \'ery small specimen \\ith arbitraiily large mag- 
netization into the originally empty stmcture. Under the heading, 
"Small Magnetization — Arbitrary Sample-Size," we shall discuss the 
propagation constant for a pencil of ferrite of any radius, coaxial mth a 
cyhndrical ^vaveguide, the space between guide wall and pencil being 
filled with an isotropic medium whose dielectric constant equals that of 
the ferrite. This is discussed in preparation for the practically more 
important case of a ferrite pencil of any radius in an air-filled guide. Here 
the unperturbed state of the system, when the pencil is umnagnetized 
and therefore isotropic, is already rather complicated and require some 
preliminary calculations. Under the heading "Small Sample-Size — Ar- 
bitrary Magnetization," we consider the case of a thin pencil of ferrite 
in an originally empty guide. 

Another topic, not easily treated except by perturbation methods, is 
that concerning end effects in samples of finite length. After a prehm- 
inary discussion of internal reflections in an extended slab of ferrite (a 
problem which can be treated rigorously), two cases are considered: a 
ferrite slug of arbitrary length, closely fitting a cylindrical guide, and a 
thin disc normal to the guide axis, of arbitrary size. In these cases in- 
terest centers around the effect of sample length on Faraday rotation, 
though for the ferrite slug a subsidiary effect, that of mode conversion, 
is also mentioned briefly. 

It should be emphasized that the perturbation methods employed here 
are not in themselves novel. They are standard to most linear eigenvalue 
problems of physics, and have been used in connection with electromag- 
netic problems by many authors. 

The remainder of the paper is devoted to a discussion of a ferrite-filled 
"cable" in plane parallel form, using the methods of Part I. The treat- 
ment is kept in terms of saturation magnetization and magnetizing field, 
and is based on Polder's equations. The paper concludes with an adden- 
dimi to Part I, which reports some calculations and graphs of field pat- 
terns in a cylindrical waveguide completely filled with ferrite. 

1. PERTURBATION METHODS 

1.1 General Method 

A number of authors have made applications of perturbation theory 
to the problems of propagation in gyromagnetic media and the exposi- 
tion which follows is included mamly for completeness. We shall develop 
the subject in the following fashion: it will be supposed that the unper- 
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turbed system is a wave guide containing a medium whose permeability 
and dielectric tensors are diagonal and isotropic, but may vary over the 
cross section of the (i:uide, although not in the 2-direction along the guide. 
For this system it will be assumed that a complete set of normal modes 
exists for which appropriate orthogonality relations are known. The 
perturbation of the system will then consist of changes in the permea- 
bility and dielectric tensors of the medium, including the addition of 
non-diagonal terms. If these changes are to be genuine perturbations, 
they must be of one of two kinds. Either, the variation in the properties 
of the medium is confined to a limited region, small in volume in some 
appropriate sense, in which case its magnitude may be large, or, we may 
have a small fractional change in the material properties extending over 
a considerable volume. The fields in the guide may be expanded in the 
noi-mal modes and a system of equations is developed for the 2-de- 
pendence of the amphtudes of these modes. These equations are then 
solved approximately, making use of the smallness of the perturbing 
terms. The results may then be specialized to the various situations of 
interest. 

Let us suppose that the unperturbed permeability and dielectric con- 
stant are ni{x, y) and ti{x, y) respectively and that the system is now 
altered so that it possesses a permeability tensor 

M2(a;, y) -Mx, y) 



jk{x, y) fX2(x, y) 




and a dielectric tensor 

fiix, y) 








M3(a;, y) 



-hi^, y) 

jv(x, y) f=2{x, y) 

sgC-c, y) 

Max%veirs equations for the perturbed system may be written, using 
the notation of Parts I and II,t in the form: 

dHt* 
dz 
dEi* 






dz 



- jwiEt — wiEt* = 0, 



(1) 



V-Hi* -jo:€,E, = 0, 
V-Et*-\- jcjfi:^, = 0. 



f We omit the vector signs from all transverse vectors, whicli are sufficiently 
labelled by the subscript "t." 
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These may be rearranged as 

dz 
V*E. - ^ + j^niHt = -Mfi2 - ni)Ht - cckHi* = Bt , (2) 

V-Ht* - JLoeiE, = ia)(e3 - ei)-E. = A, , 

V--Et* + ji^tiiH, = -jw(/i3 - Mi)-H'j = B, , 

where At , Bt , A, and 5z are introduced as abbreviations for the terras 
on the right hand side of the equations. E, and H, may be eliminated by 
substituting in the first two equations the expressions 

^^^ V-.^,* -A. 
and 

tiz = : ^ • 

The two equations so obtained are 



— V* — -\ r— + JWCi^i = -At + — V* — 



and 



iv*f^^)-^-^ + ;V.ff, = B, + lv*^. (3) 

yu \ ei / dz joi ci 

We now suppose that Et and /?[ can be expanded in the form 
-E( = X) an{z)Einix, y) 

n 

and 

Ht = T.K(.z)Hu(x,y), 

n 

where 

^.^e-^^"^ and H^ne"'^"' 

satisfy the unperturbed form of equations (3) and the boundary condi- 
tion that tangential E vanishes at the guide wall. These equations have 
solutions for certain values of /3„ only, but, if Et„ , Hm are solutions for 
^„ = c > 0, then S,„ , -F,„ are solutions for /3„ = -c. We shall assume 
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hereafter that Em and Hi„ pertain to positive /3„ values. For a given 

unperturbed mode it follows that -^^-^ reverses sign when the direction 

an{z) 

of propagation reverses. Substituting these series for Et and Ht in equa- 
tions (3), one finds 



dz 



+ jl^nttj, 






uid 



? \tI + ^'^"''"] ^"'* = "■^' - ^ "^^ 



A, 



(4a) 



(■4b) 



The orthogonality relation between functions of different n has been 



given by Adler. It is 



f E,„*-ff,„dS = 0, 



where n ^ m\ the tilde denotes the complex conjugate and the integra- 
tion goes over the cross section of the guide. We consider the fields to be 
un-normalized and write 



Clearly 



j Etn*-Ht,.dS = A„, 
f Htn*-St„dS = -An- 



(5) 



Multiplying equation (4b) by Hi„- and equation (4a) by ^i„- and inte- 
grating each expression over the cross section we have 






= -fArS,„dS-\-^ [ E,n-V*?^dS, 

Now we use the identity 

I GfS7*FdS = -f F{Grds) + f FV-G*dS. 



surface 
This yields 



boundary surface 
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J /i] Jb'dry Hi J Ml 

= jco j BM.n dS, 

since Etn-ds = 0. H^n is the H^ field appropriate to the rt**" mode. Again, 
we have 

[ Ht.-^*~= -f ^iHtn-ds}-\- f^V-S.n*dS, 

J tl Jb'dry ei •' «1 



.th 



since A, = on the boundary. E^„ is now the E^ field of the n mode. 
Making use of these results we obtain the equations for the amplitudes 
in the form: 



dz 






dS 



(6) 



Restoring the full expressions for At, Bi, A^ and B^ , we have 






I (62 - e,)ErEtn dS -jf vEt*-Ei„ dS 

+ J (m3 - ni)HM,n ds\ , 
J {m2 - t^i)HrHtn dS-jf >iH,*-Hin dS 

+ I (ea - 6i)E,^.„ rf5 



(7a) 



(7b) 



Equations (7a) and (7b) are, so far, exact, but they involve, on the right 
hand side, the functions E, , Hi , E, and H, which are still unknown. 
We are interested in those cases where the integral terms are small, 
either as a consequence of the terms (t2 - €i),v and so forth being small, 
or of then- being finite only over a small region. In the first case the fields 
El , //, , E, and H, may be replaced in the integrals by the values which 
they would have before the perturbation was made, In the second case 
this is not possible since a large change in the material constants of a 



•a- 
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region alters the field substantially within that region. Here, then, we 
hiu'e a preliminary problem to .solve, namely that of determining the field 
in the perturbed region in a zero order approximation. 

Perturbation problems may be divided into two classes by another 
distinction. The changes in material properties may be independent of 
the 2-coordinate, so that the new problem is to consider propagation in a 
unilorm guide differing slightly from the original one. Typical of such 
prol>lems is that of a waveguide containing a ferrite rod of infinite length 
parallel to the 2-axis; the perturbation consisting here of the change in 
the ])ruperties of the rod when it is magnetized. Clearly, in such cases, 
solutionsforwhich all field components vary as exp — j0z are still possible 
and the perturbation e(|uations (7) become equations to determine (3. 
On the other hand, there is a class of problems for which the perturbation 
is confined to a limited region in the 2-direction, and we are interested, 
perhaps, in the reflection and transmission coefficdents for a wave inci- 
dent upon the obstacle. Here, for example, we might think of the case 
of a disc of ferrite across the guide. If we remain in the range for \\hich 
perturbation theory is valid the changes in the amplitude of reflected 
and transmitted waves will be small, but the changes in phase may not 
be, if the perturbed region is sufficiently long. In the latter case, it would 
be possible, if the perturbation were uniform in z over the region in which 
it exists, to find solutions going as exp j^z, as described above, and to use 
these to fit the boundary conditions at the ends. It is also possible if the 
perturbed region is long, with slowly varjdng properties, to obtain suitable 
approximate solutions by the WKB method. Some of these cases \\t11 
arise in the examples which we treat below. 

1.2. Perturbations Uniform in z 

We consider first the general case in i\"hich the perturbation is uniform 
in z. In the absence of the perturbation the ?»''' mode is to be present. For 
the fields Ei and Hj. , in the perturbed region we write am(^)Etmo{x, y) 
and a„,(z)^;mo(-T'', y), respectively, where 0^(2) is ttie amplitude function 
for the »i^'' mode. If the perturbed region is one in which, for no magneti- 
zation, the material properties difTer only slightly from their unper- 
turbetl values, we may justifiably identify E,mii with Etm and 7/j„o with 
Hj„ . If the material properties are appreciably changed even in the 
absence of a magnetic field, £(mu and i^.mo have to be calculated by 
an independent method. For a„, we put ^4.„,e" , where /3 = /3,„ + 5(3 
and 5/3 is small. Similarly for Hi and E^ , we write hm Ht^o and b,„ E^mo , 
with h,„ - B^e"^^'. With such assumptions, the m"' set of equations (7) 
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gives an equation for 0, while any other set, with n ^ m, gives the 
excitation of the n^ mode. Substituting in (7) we have 



^"' ^ (8a) 



and 



701 



^- ^' (8b) 

+ (63 - eO^^^oi^z™] ds) B„ ^ jMB^. 



Ignoring squares and products of small quantities, one then has 

50 = \i (L + M). (9) 

The first example to be considered is the effect on the propagation in a 
circularly cylindrical waveguide, when a coaxial, magnetized pencil of 
fenite of very small radius is introduced. The guide radius is ro and that 
of the pencil is n ■ Before the ferrite is introduced, mi = t^o and d = eo , 
where mo and eo are the free space values. The unperturbed fields are 
those of the usual TE and TM modes in round guide. It is necessary to 
calculate first the zero order electric and magnetic fields ^Wthin the 
magnetized pencil; it will be sufficient to work out the magnetic case and 
deduce the electric one by analogy. Since the cross-section of the pencil 
is very small and transverse propagation effects consequently negligible, 
the internal field may he calculated by solving a static problem. The 
transverse magnetic field before the pencil is inserted is Htm and it is 
assumed that the pencil is so small that over a circle with a few times 
the radius of the latter, Htm is essentially uniform. We must now solve 
Laplace's equation for a cylindrical rod immersed in a magnetic field 
which is to be uniform at large distances. Within the rod, B, = fiHt 
— jicHt*, and at its surface the usual boundary conditions prevail. 
Hereafter we ^vrite n for ^2 - 

The fields are derivable from potentials *out , *in , which are of the 
form 

■l>i„ = (Ht^o-r), 

fu -\ , («-^) 
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where r = (x, y),ais a constant vector and the coordinate system has 
its origin at the centre of the rod. Continuity of tangential H at the sur- 
face of the rod requires 

HimO = Him + — r 

and then 

2 

r- 
The normal derivative at the surface of the rod is 

1 r ^* , 3$1 
f.^Bi + yeyj 

or, externally, 
. Matching the normal B's at the surface then gives 

Mo[2 Iii„, — Himo] = i-iHtmO — JkHi„,(j*, 

or 

-«imo 1—- — ^ . (10a) 

(m + Mo)'' — K- 

In a similar manner, one would find if the dielectric constant of the rod 
were e, 

E,^, = ^J^. (10b) 

The longitudinal fields E,^ and H„„ are unchanged ^Wthin the rod. 
Turning now to the expression (9) for 5/3, we have in the present-case, 

^Am Jponcil L e + Co 

(m + mo)- ~ k^ (m + md)- — k^ J 

where we have anticipated that Am is real, which we verify below. 
Since the integrand is constant the integral may be replaced by irri 
limes the value of the integrand at r = 0. We consider now a TE-mode 
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■\\ith variation e'""^, n ^ 0. We shall have, 



where 



and 



^™ - J„(»..^Je'^ J'nin.^) =0 



2 2 ''■nm 



For the fields on the axis, one finds for n = ±1, 



and 



nHr = -jH^ = -jn -^ e 



nEi = jEr = jmn ^ e' 



U \n\ ^ 1 there is no first order perturbation. We now have 



W 2 



e — eo 2 2 T^nm , 

eo —r — - w jUo — r "T 



To- M + WK + fi(i_ 



But we have 



A„. - 27r [ " Et.,*-3,„rdr, 
Jo 

j-Unm r 



Z , ./n(x)'' 



.T^ 



a; rfj;, 



= — 2t0JM0^to - " n"" ("nm " !)■ 



and then, 



5^ = 



1 ?']' ■)/n,n" 

2)3™ ro- J„(Unm)'(u„m- " 1) 



^':- 



+ n.K + Mo ti 



1 — Col 



(11) 
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For Tj\'I modes we have 



Et^ = -j/3mVx^ E, = 

Htm = iweoV*Xm 



n' 



., Xn 



where 



and 



Xm — <Jn [3nm — J 6 i <J n\Jnm) — 



o 2 2 Jnm 

1 ro^ 



Proceeding as before, we find, 

2/3m To^ Jn' ( Jnm) " L^" M + «•« + MO 



2 M + ?l« — MO I o 2 € — to 

po — r , + Mm ■ — i — 



(12) 



A problem which is of some interest, although not of immediate prac- 
tical significance, is that of a ferrite pencil of arbitrary radius and infinite 
length in a round wave guide, with the remainder of the waveguide filled 
A\'ith a non-magiietiu dielectric, whose dielectric constant, ei , is equal to 
that of the ferrite. The ferrite is supposed to be only weakly magnetized. 
For such a problem, we have, j 



S0 



^£\m •'pencil 



Htm is the field of an unperturbed TE or TM mode in the dielectric-filled 
guide, ft — tio and k are supposed small, but n , the radius of the pencil 
need no longer be small. 

For TE modes, wo have as before (again excluding the case n = 0), 



and 

nm 



^,„ = w eoMo — — ^- = Pi — 



^4-1 rnr^iriWii I ai aK ' 
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Hence, assuming the pencil coaxial with the guide. 



+ ^JjU.^]\rdr 






n Jit \ To/ 
= 2^^ J [(m - mo) £'"'^''^'" ^xJ':(x) + ^'j dx 

+ 2Kn I Jn{x)Jni-'c) dx , 

= 2jr^„," (m - Mo)f xJ^n(a:)/„'(x) +|-^?(x) 






^0 



Making use of the value of An, found in the preceeding paragraphs, we 
have 



^ 2 


LVmo j\ 

Jlix)) 


1 " r2 /^ ^A 

(ri/ro) Mo \ ro/_ 


(13) 


For TM modes. 








i^im = 


-jj3mVXm , 






Htm = 


^wei Wm , 







x„ = /.(.-.. -J.-, 



^m = to eijUo 



3tim _ o 2 

— -^ - Pi 



2 3'' 



n' 



n' 



In this case, 
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The value of A™ for this case is 

■ 2 

Am = -2lrmiP„-^ j',Kjnn.)- 



The value of 6)3 now becomes 



2 

X 



(1-i) 



We note that for a ferrite filled guide A\-ith ri = ro , the nonreciprocal part 
of 5/3 vanislies which confirms a result found in Part I of this paper for 
weak magnetization. 

The very high dielectric constant of the ferrites (about 10) puts rather 
severe restrictions on the size of the pencils to which perturbation theory 
is applicable, even for weak magnetization. This limitation would be 
substantially relaxed if we possessed exact solutions for rods of high 
dielectric constant inserted into round guide, which could be used as the 
basis for magnetic perturbation calculations. Unfortunately, the only 
extensive published calculations of this kind are for dielectric constants 
less than 3. However, at the suggestion of M. T. Weiss, a calculation of 
the propagation constant of the lowest mode varying as e*'^ in a wave 
guide contaiTiing a coaxial dielectric rod (ei = 10) has recently been made 
in the Mathematics DepartmeTit, for varying rod diameter, but for a 
single value of guide radius equal to 0.4 times the free space wavelength. 
With the aid of this information, which was made available to us, the 
magnetic perturbation calculation has been carried out. 

As before, the radius of the guide is ro and that of the rod is n . The 
dielectric constant of the rod is ci . We consider first the propagation in 
the uiimagnetized case. Since we are considering only one mode, namely, 
that with an angular variation, t;^"', and of the lowest order radially, we 
need not identify the E's and H'n by a label. We use a subscript "1" for 
fields in the dielectric and "0" for fields in the empty part of the guide. 
In general, we have 



where 



2 2 

a = ill en — 
V'E, = -a'E.. 



— -i.-_.«--»-,«.- tvSE'i^- 
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and e, m refer, for the time being, to the dielectric constant and perme- 
ability of the region considered. It ^nll be convenient to put 



// 



V Mo 



H..t; 



* - M _ -. 

- = e; — = m; 

en Mo 







a 



w V enMo 



wV€nMo 



and to measm-e lengths in terms of — 7^= . We shall continue to use 

coVMofo 

V and V* with the understanding that they refer to the scaled units. 
We now have 

_2 .- a2 

a = eti — (i , 

V'E. = -a'E.. 
At the surface of the rod E,,H^, E^ and H^ are continuous. We must 

have 



50- 



n \dr /oj ar L^i \ 






and 



where 






^^h©. + ^^^<^') 



€1 

ei = — 



ao 



= 1 - ^' and 



ri 



_ 2 - 32 

«! = ei — p . 



/i = 1, everywhere, if the unmagnetized ferrite has the permeability of 
free space. 

These equations may be rearranged in the form: 



l-iU(n)=^i. 



O-o' «i- 



ao- 



1 (dH\ ~ 



+ k 



1 (dH^ 






+ .- 



1 /dE^ 



Si^ L £^z{ri) V ar /i_ 



ii?.{ri), 



A%('-i). 



(15) 
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Withiii the fliolectric, since all fields are bounded at r = 0, both E^ and 
H^ are proportional to Ji(air) and, e\'identh', 

in the notation of Part 1. (E;)o and {H,)(, , similarly, mil be those two 
linear combinations of ./i(aor) and Yi{aor), which, respectively, vanish 
and have zero normal deri\'ative at r = ro , in order to ensure the van- 
ishing of the tangential fields there. The functions 

r d{H,), , r a(.E.)o 

and 



(H.)o Br "'^" (E.)o dr 

vn\\ be called //(ao'O and Giaor) respectively. 

Eliminating E^ri) and H^iri) from equations (15) we obtain the char- 
acteristic equation of the problem in the form 

^ ("1 _ JlY = /^Wi) _ ^(Son) \/^ F(5iri) _ Gimn) \ 

The perturbation in the present problem is that due to a mild mag- 
netization of the rod and referring again to equation (9) we have (in 
unsealed units) 

5/3 = --^ f [(m - Mf.)//r^. - >-ff,*-^r] dS, 

Avith 

An = f E*-StdS. 

Thus, in the scaled system, 



5^ = -o- 



J- Jo L\Mo / Mo 



r dr 



2 r" 



i" Et*-Hirdr 
Jo 



(16) 



The evaluation of the nonnalizing integral in the denominator is an ex- 
ceedingly tedious business and it seems advisable to avoid it. This may 
be done in the follo\™ig manner. The characteristic equation has been 
solved for numerous values of r\ and ^ may be considered to be a reliably 

known function of i\ . In particular the slope ~ is known. But we may 

dn 

also deduce -r— , by a perturbation calculation in \\hich wo start ^ith a 
an 



^^.--_. 



-J.i-A«i -- f.. 
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rod of radius n and increase the latter to n + dr^ by changing to to ei 
in the shell r, < ;■ < ;i + dr^ . For such a perturbation, since J?^ and E^ 
are continuous at the boundary, they suffer no change when the Iwundary 

moves; Er however is discontinuous and {EX = -(-Kr), where (AV), and 

to 

(Er),- are the normal fields just.outside and just inside the rod. The per- 
turbation formula (unsealed), thus gives 

&0 = :ZJ^ (e, - .o)(| E. r + I K I' + - I iE.)i r)-27rMn , 

or, scaled, with r and /'i representing scaled radii, 

op n /To ~ 

^ / Et*-H,rdr 

The formula (16) for the magnetic perturbation may now be written 

L ^ r Ut. - i) I H, r - j^ {H,*-Hd\rdT 
j.a — 7j_«!l -Jo L\Mo ' t^o __J _ 

^p - -e,- \' I EM) r + I ^.('-0 I' + ^t I i^r)- ' ' 

Inside the rod, we may write 



ft -i^ 

^' Kin) 



E, = jcE, 



and then 



The integrals are readily carried out and are as follows: 



I 



Ht \r dr 

E.'in) 



0!f 



{ii' + cY)( F{a,n) + 



F'{,dL\ri) , aiVi' — 1 



+ 



2 



- 2c^e, I , 



r {Ht*-Hi)rdr 



^ jFJin) 
Si* L 



(,^ + c^m - 2c..^(F(..rO + ^^ + ^^^')]' 
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The terai in the denominator may be evaluated by using 

r or 

dr r 
We obtain 

I EM) r + I Ee(r,) T + ii | (E.)i f 



= EM) 



2 ™.l 



The perturbation may now be written; 

*5 5n 



]■ 



il - 1 

(17) 






/ Mo 



rM^ + [^ - cF(5,ri)]= + ei[c - ^F(airi)]^ 
where 

^ = - — ■ + /' (airij + — — 

and c may be obtained from equation 15, and the definition of c 

G{mn) _ _ F(5iri) 

ao^ ^ 5i^ ai^Giaori) — eimF(airi) 



c = 



M^: - 1) 



with 

/^^_ ^ 1 _ A'"i(5oro)./i'(aori) - /i(aoro)A'"i'(aori) 
Ni{aoro)Ji{aari) — Ji{aoro)N,(aQri) 

Fig. 1 shows the propagation constant as a function of the relative 
diameter, n/ro , of the dielectric rod in the unmagnetized case, with 

€i = 10. Fig. 2 shows the derivative j} — r^ as a function of rJro . The 

rf(r]/ro) 

guide radius is 0.4 times the free space wavelength. From equation (17) 

5^ may be written as 

5^ = P(n - no) + Qk. 



--■-. f .-. t-.^i'-*?--^- rS-M-J.-.^ ,j»-rf»l,a« !<•»*»?- -.- 
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Fig. 1 — Propagation constant of a circular guide containing an unmagnetized 
coaxial, dielectric rod («/eo = 10). ri is the radius of the rod and ro the radius of 

the guide, ro = 0,4Xo , where Xo is the free space wavelength. 

The computed values of Q and P - Q are shown in Fig. 3 as a function 
of relative rod radius. P - Q is plotted since P and Q are very nearly- 
equal.* 

1.3 Perturbations Non-Uniform in z 

So far we have been concerned only with structures indefinitely ex- 
tended in the direction of wave propagation. In practice the non-recipro- 
cal element is, of course, finite, but end effects can frequently be ne- 
glected, since the element is matched at its ends (by tapering of the finite 

* 11. Seidel and Miss Rl. J. Brannon, at the suggestion of M. T. Weiss, have re- 
centlv calculated the dielectric loss for the guides containing a <lielectric rod de- 
scribed above. By combining such information with that obtained here it is possi- 
ble to discuss figures of merit (degrees of rotation loss in db) for various pencil 
radii. Such an analysis is being made by M. T. Weiss and will appear in an ar- 
ticle, by S. E. Miller, A, G. Fox and M. T. Weiss, in a forthcoming issue of the 
Journal. 
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sample, for instance). The matrhing could be accomplished in such a 
way that the transition region, whose characteristics would be very diffi- 
cult to compute, should contribute little to the overall non-reciprocal 
behavior. Therefore, in many cases, the theory for the indefinitely ex- 
tended sample is adequate. For some special purposes, however, it is 
desirable to mismatch the sample deliberately. For instance, Rowen^ 
has suggested that the change in Faraday rotation, due to internal re- 
flections in an unmatched specimen, can offset to some extent the fre- 
quency dependence of the rotation which is implied by the Polder rela- 
tions, broadening thereby the useful bandwidth of the device. 

Consider an infinite slab of ferrite, magnetized in a direction normal 
to its two parallel plane bounding surfaces. A circularly polarized wave, 
normally incident on the slab, will be partially transmitted, and, since 
for such a wave, the medium behaves as though it had an ordinary scalar 
permeability, the phase and amplitude of the transmitted portion are 
readily calculated. It is clear that, as the result of multiple internal re- 
flections, the phase of the emerging wave \y\U differ from the value ap- 
propriate to a single trip through the slab (such as would be obtained 
were the slab perfectly matched). Both amplitude and phase of the 
transmitted w^ave will depend on the electrical thickness of the slab and 
its dielectric constant, and on the effective permeability. The latter 
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P-Q 



Fig. 3 ^ Q and P — Q versus n/ro ■ 

differs for right and left circular waves, so that a plane polarized signal 
(which is the sum of equal right and left circular components) will 
emerge, in general, elliptically polarized, with the major axis of the 
ellipse tilted from the. polarization at incidence by an angle which differs 
from the single trip value as the result of internal reflections. It is clear 
that this change in rotation can be calculated in a very elementary way. 
When the sample is confined to a waveguide a similar effect occurs, 
but its calculation becomes extremely involved, at least for arbitrarily 
large magnetizations. The reason, which should be clear from Part I of 
this paper, is that the circularly polarized modes no longer have the same 
field configurations inside and outside the sample.* Therefore any inci- 
dent mode excites all of the normal modes of the ferrite; these, in turn, 
give rise to all the mode patterns of the air-filled portion of the guide. 
Even if all but one of these are cut off, the excitation modifies the phase 
and amplitude of the reflected and transmitted portions of the propagat- 
ing mode. 

* This is due to the fact tliat there is now no ordinary effective scalar perme- 
ability as for infinite geometry. 
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Thus all modes have to be included in the problem, which conse- 
quently takes the form of an infinite system of linear equations for the 
mode amplitudes. This can be solved only to some approximation whose 
genera! validity it would be hard to establish. The problem could also 
be stated as an integral equation involving a complicated Green's func- 
tion, with no greater chance of complete solution. 

We are therefore forced to restrict the problem to the ranges of mag- 
netization, or of sample size, in which perturbation theory is applicable. 
Howe\'er, we begin with a discussion of the infinite, plane, loss-free slab, 
a problem which can be solved completely, and which has some bearing 
on the perturbation problem for a slug of ferrite whose cross section 
completely fills the waveguide. 

Let the plane boundaries of the slab be normal to the 2-axis, which is 
also the dii-ection of magnetization and the propagation direction of a 
circularly polarized wave incident on the boundary z = (see Fig. 4). 
In terms of the parameters p and a- of Section 2.1, Part I, the effective 
permeability for a circularly polarized wave is 



P 



fi ± K = n = HO [ i d= 



1 ± rr 



the upper sign referring to right circular polarization. The correspond- 
ing propagation constants in the slab are then 



fi± = wV' 



CMo 



1 ± 



P 
1 ± 




Fig. 4 — NonnnUy iiiagnctizecl feiiite slab. 



^s-^tf- - ■'.Kuuir'.ri^h 



■*»i>-'"JvE(;-. 
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where t is the dielectric constant of the ferrite. If d is the thickness of the 
slab, the electrical thickness is 



where So = i^y/ti^d is the electrical thickness of the unmagnetized 
sample. Let us now confine the discussion to the right circular wave. If 
the incident electric field is taken to be e~'^°', where /5o is the free space 
propagation constant, ui -\/mo«o , the incident magnetic field \\'ill be 

— B ) 
w/xo 

and if the reflected electric field is pe^^"', the reflected magnetic field will 
be 

since /3o reverses sign. Inside the slab, the electric field consists of forward 
and backward traveUing parts Tie"'^-^' and r2e'^+^ and the corresponding 
magnetic fields are 



and 






COM+ 



Finally the transmitted electric and magnetic fields will be denoted by 
and 



^^^-y^„(.-.) 






respectively. In general p, as well as the t's, will be complex. To obtain 
T3 , we write do\\'n the equations of continuity of aU tangential fields 
across the boundaries 3 = and z = d. Since the fields are confined to a 
plane normal to the z-direction, these equations are: 

1 + p = TI + T2 , 



^ a - p) = ^ (r. - r.) 

Oifio COpo 
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and 

^ ra = ^ (r,e-'-' - r,e''-% 
These four equations in four unknowns are easily solved for T3 . Writing 



and noting that 



+ <r 



a 

= ~r . say, 
x+ 



where 



-\/b 



one finds the solution to be 



T-i = 






COS 



where 



and 



tan *+ = - ( h — ) tan Q^^ (18) 



,.-[ 



4 \x^ a 



T3 1+ = cos ^0.1;+ + 7 ( — + — 1 sin dtiXj^ 



1-1/2 



(19) 



Similar results apply to a left eii-cular wave; it is necessary only to reverse 
the signs of cr, p in the expression for a:+ . Equations (18) and (19) show 
that equal right and left circular incident waves emerge \\ith different 
amplitudes and phases; hence an incident plane polarized wave emerges 
elliptically polarized T\-ith its major axis inclined to the polarization direc- 
tion at incidence. The inchnation and the ratio of minor to major axis 



,-J- 
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are determined as follows: The right and left (iirnular fields, upon emer- 
gence, may he written in terms of rectangnlar (components (with the 
polarization of the total incident field along the x-axis) 

E--\-jE- = T-e"'', 
from which the resultant field m the a;-direction is seen to be 

E^T = ^/ + -'^V" = 1 T+ I cos (o)/ - $+) + I T_ 1 COS {ut - ^-i 
and in the y direction } (20) 

EyT = Ey^ + Ky~ = \t-\ sin (wi - $_) - | 7+ | sin {wf - <!'+)., 
The amplitude at time i, (EJ + -E„j.')''', is thus given by 
EJ -\- Eyf = lr+l'+ |r_l' 

+ 2| r+ I ■ I r_ I cos [2a - ($_ + *+)] (21) 

The major axis of the eUipse is the maximum of (E^t + Eyr ) "'ith 
respect to coi. It equals | t+ j + | t_ | and is attained at 

a = >^($_ + $4-). 

Similarly the minor axis is the minimum and equals \t+\ - | t_ |. 
The ratio of minor to major axis is therefore 

I r+ I + I r_ I • 

The angle between the .r-axis (the incident polarization direction) 
and the major axis is found by substituting c^t = H(^- + *+) i" (20)- 
This gives 



E,r = (I T+ j + I 7-_ I) cos 



$4. — '^- 



E,T = (Ir+I + lr_i) sin ^^ ^^ . 
which shows that the angle is 

~ 2 ' 

1 T I and $ are plotted versus x in Fig. 5. r+ , $+ and t- , *- at given 
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Note that x can be imaginary corresponding to cut-off in the range 
-1 < a < -1 - p (p < 0). 
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Near x = \ (corresponding either to small^ or to sufficiently large u), 
$ differs from its value ■J'o for the isotropic case by a small amount 5#. 
The rotation, to first order, is then just one half the difference between 
the 5$ for positive and negative p, <j. Writing 



:c± = 1 + 5.r± = 1 ± 



2 1 ± <r 
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and expanding equation (18), we obtain 



where ^ is defined by 



cosh * sec^ ^0 — ° sinh ■^ 

P0_ 

1 + cosh^ * tan^ do 



a=^ = .^ 



[This result holds even near 0o = (n + 3^)t where tan So = «> , as can 
be seen by expanding the reciprocal of equation (18).] The quantity 
}4^o(Sx+. — SxJ) is the rotation corresponding to a single trip through 
the sample. The actual rotation is ^(S#+ — 5$_). Hence we may define 
a rotation gain as the ratio 

g{6o , a) = j— T- — r , 

, ^ 2 ^ tan 00 - , ^ ^22) 

cosh '^ sec do " — - — smh * 



1 + cosh^ * tan^ 9o 



In many cases, flo ^ 1, that is, the thickness of the sample is much 
greater than a reduced wavelength in the specimen. Then the second term 
in the numerator of (/ is always negligible compared with the first, g then 
simplifies to 

cosh ^ 
9i 



cos^ do + cosh^ * sin^ 9o ' 



This expression is plotted in Fig. 6 as a function of 6 for various a = e*. 

For given * it has minima equal to — ; — at do = i n -{- -] ir, and max- 

cosh * \ 2/ 

ima equal to cosh ^&t6o = nwin = 0, 1, 2, • • ■). When a ^ 1 (a -^ 3 

for many ferrites), cosh * is replaced by ^e* = 3^a, and then 

1 

Ql max. ~ n^ 

2 

f/lmln = - 

a 

It is to be noted that when 6 = titt, the condition for maximum gi , 
the unperturbed reflected amplitude is zero, and the ellipticity vanishes 
to first order in 5x. 



i9-rm-ii2i,'t.\- -', ~--^-- . . _', 
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Fig. 6 — Hotation giiin, £/i(i?o), versus electrical thicloiess, do , for small mag- 
netization. 

Perturbation theory enables us to solve approximately, for weak mag- 
netization, the prolilem just considered for the case in which a right 
circular cylinder of ferrite snugly fits a cylindrical waveguide. We will 
show that, with a suitable reinterpretation of the constants, equation 
(22) for the rotation gain of the extended slab continues to apply here. 

Before magnetization, a particular right circularly polarized mode, 
say the m*'', is present in the sample. The small magnetization distorts 
the pattern of this mode and changes the propagation constant slightly. 
The distorted pattern can be expanded in the series of normal modes of 
the unperturbed material. In these expansions only the coefficient of the 
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originally present in mode will be large; all others mil at most be of the 
order of the perturbation. Denoting perturbed quantities by the superfix 
+ , we have for the distorted fields 

Htm ~ G " 2^n=\ Qmntltn y 

where pmm and g„,„ are large compared mth the other coefficients. (3^ 
and the p's and c/b are determined from equations {7a) and (7b). «„ in 
these equations is identified with e~' '" ' pmn , b„ \iith e~' " ' q„„ . Since 
all perturbation integrals involving Et and H, , E, vanish in the present 
ease, we obtain 



PnQmn PmPmn — . 



f(„ - noW^Stn dS - j j KHl*Hi. ds\ 



and 

In the firat of these equations, H^m , the perturbed magnetic field, con- 
sists of QmrnHi,,, , plus an admixture of other modes with coefficients them- 
selves of order /x — mo , «■ Therefore, to first order, it suffices to write 

in the integrand, ^\^th the result: 

PnQmn PmPnin -• mn^mm > 

ymn — Qmn , 

Pn 

where 

Jt. = ^ /" (m - n.)H,^Hu - JKHtilt.] dS. 
Elimination of pmr, gives 

The case n = m determines /3m^: 

&t: = 0l.(i - ilmi&j. (23) 

-\11 other cases give, to first order, 



jp.tlM.. . ^■,;.^.i,.,fc^. 
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which, again to first order, equals 






&l - ^l 



One of the quantities p„.™ , ?„„. , is still arbitrary. If it is required that 
the perturbed field is, to zero order in I't„ , nonnalized to the same value 
as the unperturbed field, it is readily found that one can take 

We are now in a position to consider the problem of a right circular 
mode, say the r"*, incident upon the end plane, 2 = 0, of the ferrite cylin- 
der extending to z = d. One simplifying feature of this problem is that 
the unperturbed modes inside, and the modes outside, the sample have 
the same dependence on radius since the sample fills the whole guide- 
cross-section. However the modes inside and outside may have different 
numerical coefficients. Thus, if we distinguish quantities outside the 
sample by primes, the TE,r mode can be represented outside and inside 
the sample by 

WMo 

Here sgn /3 means: sign of the propagation direction, +1 and -1 for 
forward and reverse respectively. The function ^^r is given by 

where tp is the azimuthal angle, u the guide radius, and Wr the / zero of 
JJ{x). Tor TM modes we have similarly 

E\^^ = E,^^ = Vhr , 
H,,^= -^VAV.r)sgn^., 
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where in the definition of the i^, Vr is replaced by > , the r zero of J,(x). 
Since the perturbations considered here do not couple modes ^\■ith dif- 
ferent azimuthal number, s \\ill be considered fixed hereafter, and only 
the suffix r vdW be retained. 

The field at z = —0 (just outside the ferrite) Tvill consist, not only of 
the incident field and its reflection pr , but also of other modes excited 
by the perturbation. Taking the incident E/ to have unit amplitude, we 
have, from the continuity of tangential electric fields, 

(1 + Pr)E'tr + E Pn^'i- = E (^"1 + ^^2) ^ P^^^ , 

" ' (25) 

where the t„i , t„2 are respectively the forward and backward traveling 
amplitudes of the n^^ perturbed mode. In the same way, we have for the 
tangential magnetic field: 

(1 - pr)H[r - i: pr^H't. = E (r.i - r„,)q„iHu . (26) 

The changes in sign of the coefficients of the backward waves have al- 
ready been explained in Section 1.1. Let us now suppose that the incident 
mode is the TEr mode. Multiplying hothsidcs of (25) by V*\t'r and both 
sides of (26) by V^r , and integrating over the cross-section we have, 
from the orthogonality relations of these functions, 

1 + Pr = E (^"1 + Tn2)pnr , 

j3r(l - Pr) = J^ |8r(7-„l - T,^)g'„ . 

n 

But in these series, all t except Tr vnW be small, as will be all p, q except 
p„ , Qrr ■ Hence all terms, except those for which n = r, \vi\\ be small of 
second order, and can be neglected. Further, 

Thus we obtain finally, 

1 + pr = Trl + Tr2 , 

1 f \ Prti 

1-p. = (r. -r,.)^^. 
In view of equation (23) for fir^, these equations can be ^vritten, to 



'i^:'. ^. -'■t--. 
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first order in. Itr , 

1 + pr = T;-! + Tt2 

1 - Pr = ^, \.rn - Tr2) 



where 



Mr-)- = MO 



(-1)^ 



Similarly at the output plane, z = d, the transmitted amplitude t^s 

satisfies 

Tr3 = Trie + Trafi , 

Pt,. . -. -.- (28) 



Tr3 = 



(Trie-'" - Tr^n, 



where 6 = ^r^ d. Comparison of equations (27) and (28) mth the cor- 
responding equations for the infinite slab shows that they have the 
same form; however, equations (27) and (28) hold only for small mag- 
netization. Hence only those results for the slab that relate to small 
magnetization can be generalized to the present case. One such result is 
that for the rotation gain. If we re-identify variables by 

6o^Qr = /3, d, 
we obtain, for the rotation gain, 



, . 2 „ tan 6r . , , 
cosh 4'r sec 6r — — ^ — sinh i/'r 



(29) 



"' 1 + cosh^ ^pr tan^ Or 

where 

e^' = a, = ^r/^'r . 

The formulae and conditions for maximum and minimum rotation gain 
dei-ived in the previous section, apply here also, in terms of the re-defined 
variables. 

Tlie conversion of part of the in<^ident mode to the s*^ mode in the 
transmitted wave can be examined by multiplying the matching equa- 
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tioiis by V*^« or V^, and integrating over the cross-section. It is then 
found that the transmitted ampHtude of the n^ mode is 

T = 2 Z (^^ + ^»)(cQS 9„ - cos dr) + j{l + Z„Zr){sm g„ - sin fl,) 
r„3 7r« „ (2Z„cos0„+i(Z2„+ l)sine„)(2Z,cose,+y(Z?+ l)smdr) ' 

where 

Or = ^r (i, Zr.n — —, 



fir 



and 

7rn = 



+ 






2 * 



The perturbation theory just outlined assumed a guide closely fitted 
with a slug of ferrite of finite length, and slightly magnetized. The per- 
turbation consisted in the small changes in permeabifity. Here we treat 
another kind of perturbation in which the ferrite does not fill the guide 
completely, and in which the magnetization can be arbitrarily large, 
but the sample is a thin lamuia, mounted normally to the guide axis 
which only slightlj' perturbes the field pattei'n. Its shape can then be 
considered arbitrary; its thickness will be assumed uniform and very 
much smaller than a wavelength in the material. 

Under these concUtions, the equations for the perturbed amplitudes 
of the right-circularly polaiized n^ radial mode are 

^" + JI3J>., =j^- j in - no)H,H,„ dS -j f kH,*H^. dS 

+ /(£ - eo)E,-E..ds\, (30) 
+ i)3na„ = j ^[j (e- eo)ErEt^ ds] , 



36^ 
dz 



in the usual notation. The terms on the right are assumed different from 
zero oidy in the range z,z + 5z occupied by the sample, and the integrals 
are extended over the cross section of the sample only. £",„ , Hu are the 
mode functions for the empty guide. The fields i/j , Et , E^ are the actual 
fields inside the sample. For the first-order calculation Et and Hi are 
identical mth the incident field and E^ (by continuity of the normal 

component of I*;) is - times the incident E^ . Now the incident fields may 



i^'-'tOtw-l^ -■■C' 
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be taken simply as 

(since, invacuo, &„o = a„o sgn )3„ and sgn /3„ = +1 for the incident field), 
and the amplitude a„o may be taken to be unity. Therefore equation (30) . 
may be written 

^^ = ;■ -^ r /■ (^ - Mo)//,„^,„ dS-jf KH,.*Rin dS 
dz An L"' 

+ / (l - fj eo^.n^^™ dS^ - jPnK , 
^ =jJ^ f (e - eo)E,nEt. dS - jff.a^ . 

dZ An J 

We solve these equations by integrating through the small thickness of 
the sample. Since we are interested only in results of first order in dz, 
a„ and K on the right hand side of these equations may be replaced by 
a„o , &„o ; that is by 1, 1. Writing 



^' ~ An 



dS 



/ (m - m)I^infftn dS~ 3 j KHtn*Rin 

+ J(l -f\eoE.J,„dS 
1, = -^ f (e - €o)Et.Et. dS, 

An J 

we have, upon integration from the incidence plane (1) to the exit plane 
(2), 

a„2 - Oni = j(Ia - )3„)S2, 

hm - 6nl = KIb - 0n)5Z. 

a„i , the amphtude at the entrance plane consists of a„o = 1, the incident 
amplitude, plus a small reflected amphtude p. Thus a„i = 1 + p. Simi- 
larly &„ consists of a„o sgn (3„ = 1, and the small reflected ampUtude 
p8gn^„ = —p (corresponding'to a backward wave). h„2 is equal to a„i 
(since the transmitted wave is a forward wave). Hence we have 

a.., - (1 + p) = KTa - 13.) 5z, 

a,2 - (I - P) = K^B - ^n) hz, 
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whose solutions for the transmitted and reflected amplitude are 



a.2 = l~j[^„-^-±±l£)8z 



and 



P = 3 F^ 33. 



Thus, to first order, the transmitted wave only undergoes a phase change, 
of amount 



A similar result applies to a left circular wave ; the only difference arises 
in the contribution 

to — which vnW merely change sign. The remainder of I a , and also Is 

are unchanged. Thus the Faraday rotation of a plane polarized wave, 
which is one half the difference of the phase change for right and left 
circular waves, is equal to 



ijZ-in ''oroa of diao 



The integral on the right is real, and, when the disc is circular of radius 
fi , is of exactly the type that has already been encountered in the case 
of the ferrite embedded in a material with the same dielectric constant 
(Section 1.2). Here, however, k need not be small, so long as 8z is suffi- 
ciently small. In using the results of Section 1.2 it must be remembered 
that the mode functions there related to a dielectric, while here they 
relate to air -filled guide. When this is taken into account, one obtains 
for the specific rotation of a TE-mode: 



5^ I3„ 



8Z (Wn™" - l)Ji{u„^) 

a TM-mode: 



[vA^-til (-^' 






(31b) 



^--l ..- -.^.^j.^,^,^-.5,,,i^;ji^t^j;^,^i,j^gl(;^ 
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results which show that the rotation in a thin disc is independent of the 
dielectric constant of the disc. 

2. THE PARALLEL PLANE CABLE 

Measurements on ferrites are sometimes made by a coaxial cable tech- 
nique. The fei-rite fills a section of the cable, and is magnetized axially. 
In this section we consider a simplified problem of the same kind : that of 
wave motion between parallel conducting planes bounding a ferrite slab 
magnetized m the propagation direction. The characteristic equation for 
the propagation constant has been established by a number of authors, 
particularly Van Trier,^ but as in the case of the cylindi-ical waveguide, 
no analysis of the equation in terms of the laboratory variables appears 
to have been made. We shall make such an analysis by means of the 
methods already used in Part I for the cylindrical waveguide. The dis- 
cussion will be lunited to the dominant (TEM-limit) mode, and to the 
principal features of the incipient modes (shape resonances). 

The distance between the two conducting boundaries will be denoted 
by 2a. The system of axes is as in Fig. 7; the y-ajds is perpendicular 
to the walls {located at y = ±a), the s-axis is the magnetization and 
propagation direction. The microwave fields are assumed not to vary 
along the x-axis. 

The solution of Maxwell's equation proceeds just as in the case of the 
cylindrical waveguide (Part I, Section 3), with the simplifying feature 
that one of the coordinates (x) drops out of the problem. We shall use 
the reduced notation of Part I. ^, the propagation constant will be meas- 
ured in units of ajv^> the propagation constant of the umnagnetized, 
unbounded medium, a mil denote aoi^AI^, the half-spacing measured 

in terms of the reduced wavelength -- of the unmagnetized, extended 
medium. 



CONDUCTING 
BOUNDARIES 2'a 




Fig. 7 — Parallel plane cable filled with ferrite magnetized along propagation 
direction. 



GUIDED WAVIO I'nOPAGATION THROUGH GYROMAGNETIC MEDIA. TIT 1169 

It is found that (see Part I, Section 3) the longitudinal eler^tric and 
magnetic fields are then given by 

p M^i — Ai^2 rj ■ ^1 — ^2 

^' ^ ~~A 1 — ' ^^ = J 1 r ' 

A2 — Ai A2 — Ai 

where ^'i.ziy) are solutions of 



5^1.2 



+ Xi,2Vi.2 - 0, (32) 



where Ai,2 are the roots of the quadratic (16), Part I, and where, finally, 
the x's are defined in terms of the A by equations (17a) and (17b), Part I. 
The .r-component of E is given by [see equation (22a), Part I] 



(Ai - A2)S2£. = MiP« - pl^H -{- i'hU - -^-A 



#2 



or the same expression with suffixes 1, 2 interchanged. 1 Note that the 

X component of Vip is now zero, that of V V is — ). If the boundary con- 

ditions are to be satisfied at both y = -\-a and y = —a, the solutions of 
equation (32) must be either even or odd functions of y. Since Ez must 
vanish at y = ±a, we have for the symmetric (or even) case 

, « . eos xiZ/ , X , cos X2y 

1^1 = Ai ; ^2 ^ Aa , 

cos Xll COS xafl 

and for the antisymmetric (odd) case 

fi = Ai ^2 = Ao -. — . 

sni xi*^ sui xi^ 

The characteristic etiuation for these two sets now follows from the 
fact that E^ ^ at y = zta. Expressing the A's in terms of the x's by 
means of equation (171)), Part I, and writing /SAa, 1 = Xi. 2 , we obtain 

- — -„ xifl tan xio. = -— „ X2fl tan x-n 
Aixr Aax/ 

for the symmetric case, and 

, 2 tan xi« , ' tan X''« 

AiXi ^r- = ^2X2 ^- 

Xitt X2a 

for the antisymmetric case. To bring these equatiojis into a form .suit- 
able for graphical discussion, we express the x in terms of X, (t by means 



|it.A4v-^ ^-rtw.-../ 
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of the quadratic for X, and by means of the relations between the tensor 
components of the permeability (see beginning of Section 4.11, Part I). 
We then obtain 



x-y r^' *"" "/r 



— ffXi 



- M 



(33) 



= — /i/ ^"tanai/ — 

XaT 1 - Xa^ T 1 - <rX5 

for the symmetric modes, and 

for the antisymmetric modes. When Xi , Xg pairs satisfying one of these 
equations, and the Polder relation 

Xi + Xa - o-XiX2 ^ cr + p (35) 

have been found, /3^ is given by 

/3' = -X1X2 . 

Appearances to the contrary, equations (33) and (34) describe the 
ordinary TE and TM modes when p ^ 0, regardless of a (see the dis- 
cussion opposite Fig. 3. of Part I). When a = 0, the Polder relation trans- 

forms — into ^ , so that either of equations (33) and (34) 

1 — (7X2 1 — (tXi 

can be satisfied only by 

1-Xi= 1-X2^ ^^2^ ^^ L^A'^' 



(-0' 



1 — ffXi 1 — <rX2 a- V 2/ d 

and the two alternatives respectively give 



-XiXa = /3' = 1 - n' T^ 



1\ IT 



or 

Since both equations (33) and (34) are satisfied by either the n-, or the 
(n _|_ 1,^). dependence, a problem of classification arises: "Which mode, 
TE or TM, is described, as p '-' 0, by the solutions of (33) and (34)?" 
Lea\'ing aside the question of the TEM mode we note that the following 
degeneracy exists in the isotropic case: 
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A TE mode \vith antisymmetric H. and a TM mode mth sjmimetric 
Eg have the same propagation constant 

2 / ,\2 

T / .1 



^ '^ ~hV"^y U = 0,1,2 ..•). 

A TE mode vdth symmetric H^ and a TJM mode with antisymmetric 
E, have the same propagation constant 

^' = 1-^' (n = l,2, ...). 

Since application of a very small magnetization will not destroy the 
symmetry properties, it follows that the solution of equation (34), and 
the solution of equation (33), which in the limit p = reduce to 



/3' = 1 - ; 



ihS- 



correspond to a TE-limit and a TM-limit mode respectively. Likemse, 
the solutions of equation (33) and the solution of equation (34) which 
in the limit p = reduce to 

correspond to a TE-limit and a TM-limit mode respectively. When 
p 7^ 0, the new ^ are thus given by different equations, so that the 
magnetization is seen to have removed the degeneracy of the isotropic 
case. 

In the special case of the TEM limit mode (/5^ = 1 when p = 0) only 
one of equations (33) and (34) has real roots. For, when p = 0, the 
Polder relation gives 

1 - Xl' ^ 1 - \2^ 

1 - ffX, 1-0X2' (36) 

= 0, 

if )3^ is to equal unity. But equation (36) satisfies equation (34) only 
[equation (33) would demand Xi = X2 which is impossible for real $\. 
Thus the TEM mode exists as the limit of an antisymmetric mode only. 

In fact it is easily shown that for a value of a too small f less than -- j to 

admit any except the TEM mode in the isotropic medium, the only solu- 
tions of equation (33) for general a, p describe incipient modes. 



■*"-*''--■ -^-■■•*«i*■■a--■K-^^«*•--*i•i*iB■■ ..t*-» ^-J.-\i 
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Thus for rt < ^ the analysis may be confined to equation (34) which 

mil give the coiu-so of the TEH limit mode. The graphical analysis 
follows closely tliat of the ryliiidrical waveguide (Part I, Section 4.1); 
the part played by the G-function there is now played by the function 



yf^t='°«'/F^ 



L = \a/ ^ ^tana/i/ ^ 



A contour map is drawn of the function L in the upper half of the \-a 
plane. Only the upper half plane need be considered smce the L-equation 
(34) is unchanged by reflection in the origin. This, incidentally, means 
that the present .situation is reciprocal, so that cr, p hereafter can be 
considered positive. Pairs of values Xi(>0) and X2(<0) satisfying the 
L-equation can then unmediately be read off, but while they will give 
^^ = X1X2 , they do not necessarily, for given o-, p, satisfy the Polder re- 
lation, equation (35). To ensure this, the quadrant, X < 0, (t > 0, 
and the L-contours therein are transformed on to X > 0, <7 > by the 
Polder relation for fixed p: 

X. . ^ +P-^^ = T(XO. 
1 — aXi 

The .'surfaces Li = L(Xi , a) and La = L(T(\i), a) will intersect in a 
number of curves, along whose base curves in the X-o- plane both the 
Polder equation and the L equation are satisfied, and along which (3^ = 
-X1X2 is known as a function of a, p. The zero and infinity curves of 
L(X, a) are denoted byO, / when X > 0, and byO', I' whenX < 0. The 
transforms of 0', /' onto X > are denoted by (0')r , {I')t. The suffix 
n denotes the infinity (zero) eur\^es corresponding to 

1 _X= / , iV/ /I -X^ ,^'\ _ , 

= " + H N^ ; -. ^ = " ^ " mtegral. 



1 - ffX V 2/ a^ \1 - ffX 

The lines X = 0, +1, — 1, are all zero curves denoted by 0^ , Ob , Ob', 
respectively. The fine Xa = 1 is a conditional infinity curve, called /,- . 
It is an / curve when viewed from (tX > 1 for X < 1, and when viewed 
from ffX < 1 for X > 1 . Otherwise (for ctX < 1 , X < 1 and ffX > 1, X > 1) 
it is a limit curve of all possible cur\'es L = const., where the constant 
takes on any value indefinitely many times. (See Part I, Section 4, 
where the curve ffX = 1 is a conditional zero <'urve, 0, .) Fig. 8, drawn 
for a '-^ 1, p '^ 0.5, shows the part of the first (|uadrant allowed by the 
Polder relation divided into regions of like and unhke sign of L(K, <j) 
and L(T(Xi , a)) by the various 0, / and (0')r , {I')t curves. The un- 
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shaded regions are regions of like sign, and all these carry solution 
curves (dotted lines), by the same reasoning as was employed in Part I, 
Section 4.11. Two branches of the TEM limit mode exist; in the area 
bounded by (O.i) r , Og , {Ob') t and in the area Ob, /i , (Ob') t . The branch 
in the first region begins at o- = \\ith Xi.Xa given by 

AiVl - Xi^ tan aVl - Xi^ = XaVl - M^ tan d\/l - \^% 

Xi + X2 = p 

and ends at the intersection of (O^)t., Ob with a = 1 — p and jS^ = 
(since X2 = on (Oj) r). The branch in the region Ob , Ii , (Ob') t begins 
with 02 = CO , ff = 1 and proceeds towards o- = oo , |S^ = 1 . 

The region bounded by 7i , (0^) r , /^ contains an infinity of solution 
curves, the incipient modes. The n"" of these begins at ff ^ 1, Xi = 1, 
the intersection of /„ and h (the line Ic is also the transform of X2 = — 00 , 



(O'b)t 




Fig. S — Znro ;iih1 infinity i-ontoui's of L(X, , a) and L(T(\,),rT) in the first 
quiiiininl <»f llio X - it ])ljinp. Cross-hatched regions are excluded by the Polder 
reliitioji; sliiided ureas ore regions of unlike sign of the two functions. Doited 
curves are solution curves. 



'■-',T\ ■.--..., -^-.-«*ji^;-»: 
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along which Li = oo); crosses from one allowed region to the next 
through the intersection of 0„ with (Os') r and proceeds to its cut-off 
point (0„ - (O^)^). 

Formulae giving /3^ in the neighborhood of special points on the solu- 
tion curves, or for special values of the parameters, are obtained exactly 
as in Part I, by expansion of the L-equation. The results are tabulated 
below. They relate to antisymmetric modes only. 

Cut-off of the TEM-limit mode {a = I - v, &^ = ^) 



Near cut-off 



da / _ tan d\ (37) 



/ , tan d\ 



For p > 1, the TEM mode has no lower branch. 
Behavior near resonance (a- = 1): 



2 
IT 



TEM mode: /3' = -r-^p, 

n^^ incipient mode: /3^ = \ ^ — 1 h _ . 

(n = 1,2, •••) 
Cut-off of the n''' incipient mode (parametric representation) : 

p = ^ - X,; ^' = (39) 

dl ^ r^/n^ _ iVoth^j/rf^' - l),-+ 2tan_%- - el 

(The reader may note the similarities of these formulas to those 
for the cylindrical waveguide). 

Spot-point for the n^ incipient mode: 

^=K^-.i)A]+.4^' (40) 

p = (I 1- £r)(^' - 11. 
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Approximate formulae for all antisymmetric modes at small p, 
TM modes: /3' = )3„,i' v^^ , (41) 

1 — 0"' 

TE modes: ^' = ^n.'U - ^^X ^^2) 

where 



2 2 

/?..! 1 ^, 

1 


71 = L2. 


and 


n = 0, 1, 2, 


PnA i ., 


TEM mode: 




1 — (7" 




TEM mode for large a: 




^'-+-:+«&> 




TEM mode for very small a: 





(43) 



(44) 



for sufficiently small a, the /3 for the TEM mode is independent 
of a. For general a, p, except a = 1, 

f = ;-^-+^r . (45) 

1 — o-p — 0- 

Reference to the expression for n and k in terms of a, p wHl sho^^' that 
the wave progresses as though the medium had a scalar permeability 
ii^ti given by 

1/1^1 



fi,.i{ 2\ti + K n — kJ' 

3. FIELD PATTERNS IN THE FILLED CIRCUIiAR GUIDE 

In Part I we discussed propagation in a circular guide filled ^\^th 
fcrrite and longitudinally magnetized. Formulae or the field com- 
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ponents were given there, and it is of some interest to examine numeri- 
cally or graphically the distribution of E and H. Since all field compo- 
nents vary {in complex notation) as e''-"^-^"'\ with n = 1 in the present 
case, it is clear that the field patterns are stationary in a system rotating 
with angular velocity, — oj. Writing $ for the angular varialile, 6 + oit, 
in this system, one finds the following formulas: 

Er = (E+ - EJ) sin $, 



where 



E^ = -E- sin 2*, 

A'„ = E_ cos 2^ -\- E+, 

E: = Ea cos $, 

H^ = H- cos 2$ + i/+ , 

Hy = 11- sin 2*, 

//z = Hq sin *, 



E-{t) = 



Eg = (E+ + ^_) cos *, 

Hr = (H+ + HJ) cos *, 
He - (H^ - //+) sin *, 



(46) 



2xi^i(xi^o) 
1 



E+ir) = 



Mi- 



2xiJ"i(xi^o) 



J^ixir) - 



Joixir) - 



1 + 



H-ir) = - 



1 - Xi 



2xiJi{xira) 



Mxir) - 



(47) 



-1 + 



HJr) = 



1-Xi^ 

XI 



2xi/i(xi^o) 



/oCxi'O - 



//,(,.) = _i -/^(^.O 



[V 



Xi Jiixii'o) 

The terms in square brackets are in each case the same as the cor- 
responding unbracketed terms, Xi and X2 replacing Xi and xi ■ The quan- 
tities E+ and E^ are the amplitudes of the left-handed and right-handed 
components of ch'cular polarization into which the transverse ^'-field 
may be resolved at each point. //+ and //_ have a similar significance 
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for the trans\'erse //-field. This may be readily verified by examining 
the vectors (E^ + jEy)e~'"' and {/^i + jHyje^'"', which represent the 
transverse field vectors in the laboratory system. The transverse fields 
are elliptieally polarized at any point and the ratio of minor to major 
axes of the elhpses are 

\\E^\ -\E-\\ , II i/+ I - I H_ 1 I 

and 



\\E+\-\- \E^\\ 1 I //+ I + I //_ I I 

The fields so far are normalized only by the choice of a simple form for 
the function, Eo{r) ; all components may, of course, be multiplied by the 
same arbitrary constant. There is some virtue to a normalization based 
upon power flow. The power flow is given in unreduced units by 

[ {EX H% dS. 

Using the scaled units of part T with 



and // replaced by A/ — // to give it the same dimensions as E (e is 

the dielectric constant of the ferrite), the power flow becomes in the 
present variables 

P = - a/- -^ f (AV^+ + E.H.) r dr, 

V Mo WMof '' 

We shall normalize the E aiid // fields by dividing the values given by 
e(iuations (47) l)y /"'. This makes the power per (isotropic wavelength 
in the ferrite)' a constant. In Fig. 9 we show the normalized fields for 
the TEii-limit and TMii-Hmit modes as a function of ?■ for the case 
fo = 5.75, \ p\ ^ 0.(> and Ke\'eral values of a. The behavior of the modes 
as a function of a and /) for this radius is shown in Fig. 14 of Part I. 
We also show the amplitudes for the isotropic cases, a = p = 0. It may 
be recalled from the discussion of cut-off points in Part I that, for the 
TM mode at this railius, when cut-off is reached at a = —0.4, the am- 
plitude of the H^ field is o\erwhelmingly greater than that of the others. 
E\en when normalized to the same power flow, the field amplitudes 
for a given a are undetermined to a factor of ±1. This factor has been 



'^lodk^ I'c .^a^ --V^-.- -.-a- •■ ■ ...J.„ ,, 



0.4 





-\ 


^""^•^ 
















TE„ 












--^ 










p = -0.6 














"^^ 


^ 


























X 


N.. 


























E*^ 


\ 








S 0.1 

13 
















M" 




^ 





— 










^.^-- 


.— ■ 





E° 






\ 


^.,^ 


a. 


^,_ 


_,-s; 


ti^ 


' 


"^ 






H' 


^ 


ir^ 


^^^^'^ 
















E" 


- 


— ■ 


..* 


,11''' 




















^^^::L 


















>- 


,"' 




1 




-a2 












if 


.^■^ 
















„, 


-— — • 


•■** 
















-a3 
























H-. 




2.5 3.0 3.5 

RADIUS 



Fig. 9 — Amplitude of the field components in lillod cylindrical guide, longi- 
tudinally magnetized, as a function of radius, for tlie TEn-liniit and TMn-iimit 
modes at several values of a-, ro = 5.75. | p | = 0-6, except for Figs. 9(f) and (p) 
which arc the isotropic cases, a ^ p = 0. The fields are normalized to a constant 
power flow per (free space wavelength)^, (a), top, (b), bottom. 
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Fig. 9(c) — See Fig. 9. 
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Fig. 9(h), top, and (o), bottom — See Fig. 9. 
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Fig. 9{q) —See Fig. 9. 
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Fig. 9(s) —See Fig. a. 
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chosen for each pattern according to the following considerations. For 
large | a \ the amplitudes should tend to those of the isotropic case ; for 
I o- 1 < (To the amplitudes should develop ^\ith increasing [ a \ away from 
the isotropic amplitudes. Thus, if the phase of the isotropic amplitudes 
is first fixed, that of all the others may be fixed by the conditions that 
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Fig. 10 — TmnHverse field patterns in a rotating aystem for four particular 
cases of Fig. 9. For any pair of E and H patterns the length of the arrows is propor- 
tional to (ield strength, but the patterns for dirt'ercnt ir-valuea have not been 
given a common normalization, (a), top, TEn-limit mode, a = 0.6, p = 0.6; (b), 
bottom, TE,i-Iimit mode, <r = -0.6, p = -U.6; (c), TM„-limit mode, a = 0.6, 
p = 0.6; and (d), ff = -0.3, p 0.6. 



GUIDED WAVE PROPAGATION THROUGH GYROMAGNETIC MEDIA. HI 1193 

they shall fall into an orderly sequence. In Fig. (10) the actual transverse 
E' and H- field patterns are sho^^'n for some representative cases. 

It will be noted from Fig. 9 that for both TE- and TM-limit modes, 
the amplitude patterns for < a < ao resemble those for a- < — 1, 
while those for — tro < o" < resemble those for 1 < a. Further, the 
patterns in the first two ranges of a are quite similar to the isotropic 
patterns; those of the latter ranges depart markedly from the isotropic. 
All of these similarities are more clearly seen for the TE modes than for 
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Fig. 10(c), top, and (d), bottom — See Fig. 10. 
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the TM modes. In Part I where the propagation in guides of very large 
radius was briefly examined, it may be recalled that in both of the ranges 
< a < cq and o- < — 1 the field patterns approached those of a plane 
wave rotating in the same sense as the pattern as a whole. Again, it was 
found that for — o-q < cr < or o- > 1, the field at most points in the 
guide was locally rotating in the opposite direction to that of the whole 
pattern. The points of similarity to the present case are clear and it is 
also evident that the TE mode more nearly approaches t^e large guide 
situation because ro = 5.75 is much further above the cut-off radius for 
this mode than it is beyond the cut-off for the TM-mode. 
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